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We  present  quantum  scattering  calculations 
for  the  collisional  relaxation  rate  coeffi¬ 
cient  of  spin-polarized  87Rb (f=  2,m  =  2) 
atoms,  which  determines  the  loss  rate  of 
cold  Rb  atoms  from  a  magnetic  trap.  Un¬ 
like  the  lighter  alkali  atoms,  spin-polar¬ 
ized  87Rb  atoms  can  undergo  dipolar  relax¬ 
ation  due  to  both  the  normal  spin- spin 
dipole  interaction  and  a  second-order  spin- 
orbit  interaction  with  distant  electronic 
states  of  the  dimer.  We  present  ab  initio 
calculations  for  the  second-order  spin-or¬ 
bit  terms  for  both  Rb2  and  CS2.  The  correc¬ 
tions  lead  to  a  reduction  in  the  relaxation 
rate  for  87Rb.  Our  primary  concern  is  to  an¬ 
alyze  the  sensitivity  of  the  87Rb  trap  loss 
to  the  uncertainties  in  the  ground  state 
molecular  potentials.  Since  the  scattering 
length  for  the  3X+U  state  is  already  known, 
the  major  uncertainties  are  associated 


with  the  X'S+g  potential.  After  testing  the 
effect  of  systematically  modifying  the 
short-range  form  of  the  molecular  poten¬ 
tials  over  a  reasonable  range,  and  intro¬ 
ducing  our  best  estimate  of  the  second-or¬ 
der  spin-orbit  interaction,  we  estimate 
that  in  the  low  temperature  limit  the  rate 
coefficient  for  loss  of  Rb  atoms  from  the 
/=  2,m  =  2  state  is  between  0.4  X  10-15 
cm3/s  and  2.4  X  10“ 15  cm3/s  (where  this 
number  counts  two  atoms  lost  per  colli¬ 
sion).  In  a  pure  condensate  the  rate  coef¬ 
ficient  would  be  reduced  by  1/2. 
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1.  Introduction 

The  recent  observation  of  Bose-Einstein  condensation 
(BEC)  in  magnetically  trapped  alkali  atoms  [1,2,3]  has 
brought  to  completion  a  15  year  attempt  to  achieve  BEC 
in  a  weakly  interacting  atomic  system  [4].  The  success 
of  BEC  in  both  S7Rb  [1]  and  23Na  [3]  and  evidence  for 
BEC  in  7Li  [2]  were  remarkable  achievements  brought 
about  by  the  development  of  laser  cooling  during  the 
past  decade,  the  design  of  optical  and  magnetic  traps  for 
holding  cold  atomic  samples,  and  most  recently  the  de¬ 
velopment  of  evaporative  cooling  techniques  to  cool 
atoms  below  the  recoil  limit.  This  experimental  success 
has  renewed  the  interest  in  collisional  loss  rates  for  spin 
aligned  alkali  systems,  since  the  binary  and  ternary 
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collisional  rates  will  limit  the  lifetimes  of  the  experi¬ 
mental  condensates.  Experimentally,  h7Rb  is  trapped  in 
the  (/  =  2,  md  =  2)  state,  designated  the  doubly  polarized 
state  or  the  stretched  state,  for  which  binary  hyperfine 
changing  collisions  could  be  the  dominant  loss  process. 
Here  /  and  ma  designate  the  quantum  numbers  for  the 
total  angular  momentum  of  the  S7Rb  atom  and  its  projec¬ 
tion  on  some  convenient  space-fixed  axis.  On  the  other 
hand,  23Na  atoms  have  so  far  been  trapped  in  the  (/  =  1 , 
ni;,  =  —  1)  state,  which  is  theoretically  expected  to  be 
more  resistant  to  binary  collisional  loss.  Hence,  the  life¬ 
time  of  the  23Na  condensate  will  probably  be  limited  by 
three-body  rates,  although  a  condensate  of  stretched 
state  23Na  atoms  may  be  affected  by  binary  collisional 
loss  as  well. 

The  purpose  of  this  paper  is  to  provide  the  most 
accurate  calculations  possible  of  the  binary  collision 
rates  for  all  inelastic  hyperfine  scattering  processes 
which  can  contribute  to  the  loss  of  spin-polarized 
ground  state  87Rb  atoms  at  temperatures  (T  <  1  p,K) 
associated  with  the  recent  experimental  observation  [1] 
of  Bose-Einstein  Condensation  (BEC).  This  spin-relax¬ 
ation  is  due  to  the  following  processes: 

87Rb(5s  ,/a  =  2,  ma  =  2)  +  87Rb(5  s  ,/  =  2,  mb  =  2) 

-4  87Rb(5s,/a,ma)  +  87Rb(5s,/b',mb').  (1) 

Having  all  87Rb  atoms  in  the  stretched  state  with  ma  =/a 
is  ideal  for  Bose-Einstein  condensation,  since  inelastic 
collisions  between  such  stretched  states  have  very  small 
rate  coefficients.  The  entrance  channel  in  Eq.  (1)  has  a 
spin/=/a  +  /b  of  magnitude /=  4.  Since  the  87Rb  atom 
has  nuclear  spin  3/2,  this  entrance  channel  can  only 
project  onto  the  triplet  a3S+u  state  of  the  atom  pair 
which  has  total  electron  spin  5=1.  Consequently  a  sim¬ 
ple  spin-exchange  model  of  the  collision  [5,6]  shows 
that  stretched  states  do  not  relax  during  the  collision. 
Our  concern  is  determining  the  small  but  significant 
rate  coefficient  for  the  trap  loss  processes  indicated  in 
Eq.  (1)  that  occur  when  the  degeneracy  of  the  molecule- 
fixed  projection  1/21  =  0  or  1  of  the  5  =  1  triplet  potential 
is  broken  by  relativistic  forces.  This  splitting  leads  to 
spin-relaxation.  We  use  standard  quantum  scattering 
methods  to  calculate  the  spin-relaxation  event  rate  coef¬ 
ficient  /Cvem,  defined  by  Stoof  et  al.  [7],  summed  over  all 
/,ma  ,/b',fflb'  channels  that  lead  to  loss  of  trapped  atoms, 
namely  channels  for  which  either  /'  and/or  /'  #  2.  Such 
collisions  lead  to  loss  of  both  atoms  from  the  trap  be¬ 
cause  of  the  large  kinetic  energy  release  (equal  to  one  or 
two  units  of  ground  state  hyperfine  splitting  shared 
equally  between  the  atoms).  The  total  rate  of  spin  relax¬ 
ation  in  a  Maxwellian  gas  (number  of  atoms  lost  per  unit 


volume  per  unit  time)  is  —2KeveMn2,  where  n  is  the 
density  of  /  =  2,  ni  =  2  atoms,  since  two  atoms  are  lost 
per  event  represented  by  Eq.  (l)._If  a  condensate  is 
present  //vent  is  multiplied  by  (2  —  £2)/2,  where  £  is  the 
condensate  fraction  [8], 

Two  atomic  parameters  are  required  to  describe  the 
separated  atoms:  the  isotopic  mass  m87  =  158425. 8we 
(where  me  is  the  electron  mass  =  9.109  X  10-31  kg)  and 
the  6834.683  MHz  splitting  between  the/  =  2  and /  =  1 
hyperfine  components.  Assuming  that  the  molecular  hy¬ 
perfine  Hamiltonian  can  be  adequately  represented  by  a 
unitary  frame  transformation  of  the  asymptotic  atomic 
hyperfine  Hamiltonians,  the  accuracy  of  the  loss  rate  is 
basically  limited  by  the  accuracy  of  the  molecular  inter¬ 
actions  we  incorporate  in  our  close-coupled  scattering 
codes.  To  perform  the  dynamic  calculations,  we  require 
four  accurate  molecular  potentials  Vs,«(R):  one  defined 
by  the  ground  A  //  state  with  5  =  0  and  with  molecule- 
fixed  spin  projection  //  =  0;  and  three  defined  by  the 
lowest  a3X+u  state  with  5=1  and  spin  projections 
<1  =  0,  ±  1 .  These  potentials  take  the  following  asymp¬ 
totic  form  [9,10,11], 

Vo,o (R)  ~  -  Cexce^  -  ( CeR ~6  +  CsR's  +  C 10/T10)  (2a) 

Vi.„C R)  ~  +  Cexce"^  -  ( CeR~b  +  CsR "8  +  Ci0R~ 10) 

+  Vssn=0(R)  +  Vson=o(R)  (2b) 

Vi.  ±  1  (R)  ~  +  Cexce-^  -  (CeR~6  +  c8/r8  +  C10/T10) 

+VsV±i(/0+HsV±i(R)  (2c) 

where  Vssn=0(R)  =  ct2R~3  and  Essfl=  ±  ,(/?)  =  -  1/2 ct2R~3 
(a  ~  1/137  is  the  fine  structure  constant)  are  the  famil¬ 
iar  spin-spin  dipole  terms  that  are  primarily  responsible 
for  dipolar  relaxation  in  hydrogen  and  the  lighter  alkalis. 
These  spin-spin  dipole  terms  have  been  elegantly  treated 
in  a  series  of  papers  by  Verhaar  and  collaborators  in 
Eindhoven  [7,12],  The  second-order  spin-orbit  terms 
Vson  which  we  have  included  are  less  well  known  in 
atomic  collision  physics,  although  they  are  well  recog¬ 
nized  as  significant  terms  in  molecular  spectroscopy 
[13,14].  These  terms,  which  are  induced  by  spin-orbit 
interactions  mediated  through  distant  electronic  states, 
mimic  the  effect  of  the  direct  spin-spin  terms  by  intro¬ 
ducing  a  splitting  in  the  {}  =  0,  ±  1  projections  of  the 
5=1  state  and  can  significantly  modify  the  spin-relax¬ 
ation  rates.  These  terms  will  be  discussed  in  detail  in 
Sec.  5,  where  we  will  show  that  they  are  of  opposite  sign 
to  the  u-R  3  terms  and  tend  to  diminish  the  spin  relax¬ 
ation  rate  for  87Rb. 
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Our  goal  in  this  paper  is  to  systematically  assess  the 
uncertainties  in  the  spin-relaxation  rates  that  are  intro¬ 
duced  by  various  uncertainties  in  the  molecular  parame¬ 
ters  that  enter  into  Eq.  (2),  and  provide  realistic  bounds 
on  the  possible  range  of  the  loss  rate  coefficient  Ktve 
that  can  be  expected  for  Rb2.  All  our  rates  are  calculated 
using  molecular  potentials  obtained  from  ab  initio 
MCSCF  codes  employing  the  highly  accurate  ab  initio 
pseudopotentials  of  Krauss  and  Stevens  [10]  for  R  <  20 
a0  (1  a0  =  0.0529177  nm),  and  joined  on  to  the  well-an¬ 
alyzed  long  range  dispersion  potentials  for  the  diatoms 
[11]  for  R  >  20  aa.  Because  of  the  usual  limitations  of 
pseudopotentials  and  the  typical  convergence  properties 
of  ab  initio  calculations  combined  with  the  extraordi¬ 
nary  demands  we  make  on  the  required  accuracy  of  the 
molecular  potentials  at  ultra-cold  collision  energies, 
these  potentials  can  only  serve  as  excellent  initial  esti¬ 
mates  of  the  short  range  portions  of  the  Vs ,n(R)  poten¬ 
tials  in  Eq.  (2). 

Collaborative  work  combining  photoassociative  spec¬ 
troscopic  data  from  the  Texas  group  with  theoretical 
analysis  by  the  Eindhoven  group  [15,16]  has  done  an 
excellent  job  of  characterizing  the  a  3XU+  potential, 
which  controls  the  entrance  channel  dynamics  for  the 
spin  relaxation  described  by  Eq.  (1).  In  Sec.  3,  we  exam¬ 
ine  the  sensitivity  of  the  scattering  length  to  variations 
in  this  potential.  We  will  accept  the  analysis  of  the  Texas/ 
Eindhoven  group  and  have  insured  that  our  potential 
reproduces  both  the  scattering  length  Ai  ~  +  1 10  a0  for 
the  87Rb  isotope  and  A ,  ~  —  300  a„  for  the  85Rb  isotope. 
In  addition,  we  introduce  a  useful  new  way  to  associate 
the  scattering  length  with  the  binding  energy  of  the  last 
bound  state  in  an  attractive  molecular  potential. 

In  Sec.  4,  we  demonstrate  a  strong  sensitivity  of  the 
relaxation  rate  to  the  shape  of  the  A  'X/  singlet  poten¬ 
tial,  when  the  a  3X„+  scattering  length  is  kept  fixed  at  its 
known  value.  In  this  case  the  sensitivity  is  not  due  to  the 
a  3XU+  entrance  channel,  but  depends  on  a  final-state 
close  coupling  effect  in  the  exit  channels.  The  sensitiv¬ 
ity  to  the  potential  leads  to  a  factor  of  six  uncertainty  in 
the  rate  coefficient,  which  is  analyzed  using  generalized 
MQDT  theory  [17]  and  especially  the  associated  half 
collision  amplitude  version  [18]  of  the  theory.  We  also 
describe  the  interplay  between  the  spin-spin  (SS)  and 
second-order  spin-orbit  (SO)  contributions  to  the  relax¬ 
ation  rate. 

We  previously  had  speculated  that  the  SO  terms 
would  modify  Rb  spin-relaxation  [19].  In  Sec.  5  we 
present  new  ab  initio  calculations  of  these  SO  terms  for 
Rb2  and  Cs2.  We  also  give  our  current  understanding  of 
the  uncertainty  range  of  the  spin-relaxation  rate  of  Rb. 
Finally,  a  summary  of  our  results  is  presented  in  Sec.  6. 


2.  Scattering  Theory  of  Ground  State 
Alkali  Atoms 

In  a  field-free  collision  both  the  total  angular  momen¬ 
tum  F  =/,  +/i,  +  (=f+  f  and  the  total  parity  p  =  ±  1 
are  constants  of  the  motion.  The  parity  is  the  symmetry 
associated  with  the  inversion  of  all  electron  and  nuclear 
space-fixed  coordinates  through  the  center-of-mass  of 
the  dimer.  The  collision  loss  rate  coefficient  Kevmt(E)  for 
total  collision  energy  E  is  expressed  in  terms  of  sums 
over  F  and p  involving  the  scattering  matrices  S (F ,p ,E) 
[8],  Each  S  -matrix  is  derived  from  a  multichannel  wave- 
function  calculated  from  standard  close-coupled  codes 
for  a  given  F ,p ,  and  E,  using  an  expansion  in  a  channel 
state  basis  \F,M,p  ),  which  describes  the 

asymptotic  properties  of  the  separated  atoms, 

(EM)  =  ,  X  I F,M,p-,y'  >  FUE,R),  (3) 

y—EJj  aJ  b 

where  f  is  the  angular  momentum  (partial  wave)  quan¬ 
tum  number  of  the  interatomic  coordinate  R ,  /  repre¬ 
sents  the  magnitude  of  the  channel  angular  momentum 
/=/a  +  /b,  and  y  gives  the  spin  channel  in  which  the 
collision  starts.  The  +  indicates  normal  scattering 
boundary  conditions  for  an  incoming  state  in  channel  y 
and  outgoing  spherical  waves  in  channels  y'.  The  chan¬ 
nel  states  are  symmetrized  with  respect  to  interchange 
of  the  identical  nuclei.  One  consequence  of  this  sym- 
metrization  is  that  odd  partial  waves  are  missing  from 
Eq.  (3)  for  a  collision  of  two/,  =  2,  =  2  87Rb  atoms. 

In  the  absence  of  spin-spin  and  second-order  spin-or¬ 
bit  interactions,  the  molecular  Hamiltonian  of  two  col¬ 
liding  ground  state  alkali  atoms  in  (ns)  orbitals,  which 
can  be  expected  to  have  zero  total  electronic  orbital 
angular  momentum  L  =  (/,  +  /,)  =  0,  possesses  two  ad¬ 
ditional  almost  good  quantum  numbers.  These  are  £  and 
/.  Although  this  is  not  true  when  one  or  more  of  the 
atoms  possesses  electronic  orbital  angular  momentum, 
such  as  an  alkali  in  its  first  excited  (up)  orbit,  it  is  an 
excellent  approximation  for  the  two  colliding  Rb(5s) 
atoms  in  Eq.  (1).  The  physical  reason  is  that  for  L  =  0 
there  are  no  electrostatic  interactions  that  cause  locking 
of  the  electron  spin  angular  momentum  of  the  system  to 
the  internuclear  axis.  Ultimately,  weak  spin-spin  dipole 
(SS)  and  second-order  spin-orbit  (SO)  interactions  cause 
the  total  electronic  spin  S  =  (sa  +  ,vh)  to  couple  to  the  axis 
and  lead  to  the  small  energetic  splittings  between  the 
molecule-fixed  spin  projections  fl  represented  in  Eq. 
(2).  For  the  moment,  if  we  ignore  these  latter  interac¬ 
tions  the  fl  projections  are  perfectly  degenerate  and  we 
can  easily  transform  the  asymptotic  channel  states 
\F  ,M  )  into  a  basis  defined  by  the  total  elec¬ 

tron  spin  angular  momentum  S  (S  =  sd  +  sb)  and  the  total 
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nuclear  spin  1(1  -  ta  +  ib),  where  ,va  and  ,vb  are  the  atomic 
electron  spin  angular  momenta  and  ia  and  ib  are  the 
nuclear  spin  angular  momenta  (for  'Rb,  ,ya  =  ,vb  =  1/2 
and  ia  =  ib  =  3/2).  This  transformation  is 

\F,M,p;-e,f,S,I)  «  2  V (25+  1)(2/+  l)(2/a+  l)(2/b+ 1) 

/a/b 

V  i  -d  ,  /a  1 

«b  *b  /b  1  I T,M,p;/,/,/a,/b  )■ 

S  /  /J 

In  this  new  representation  two  adiabatic  Born-Oppen- 
heimer  potentials,  uniquely  identified  by  the  quantum 
numbers  5  =  0  (i.e.,  the  X 1  X+g  state)  and  5=1  (i.e.,  the 
a  3S+U  state),  appear  on  the  diagonal  of  the  Hamiltonian 
matrix.  In  the  absence  of  SS  and  SO  interactions  the 
block  of  5  channel  states  and  the  block  of  I  channel 
states  are  diagonal  and  only  couplings  involving  a  simul¬ 
taneous  change  in  I  and  5  are  introduced  by  the  hyper- 
fine  interactions.  These  coupling  are  constrained  to  sub¬ 
blocks  which  insure  that  f  =  I  +  S  is  conserved,  and  we 
find  both /and  /  remain  perfectly  good  quantum  num¬ 
bers  at  all  distances.  Of  course,  at  small  distances  the 
exchange  splitting  between  the  molecular  potentials  is 
large  compared  to  the  hyperfine  splittings  and  hyperfine 
coupling  is  negligible.  However,  as  we  shall  see,  at  dis¬ 
tances  of  the  order  R  ~  20  a0  to  40  aD,  the  hyperfine 
interaction  becomes  important  and  the  I  S  coupling 
drives  the  system  back  into  the  asymptotically  diagonal 
basis  of  channel  states  I F,M,/?;/,/,/a,/b  ). 

As  seen  in  Eq.  (2),  both  SS  and  SO  interactions  pro¬ 
duce  an  energy  splitting  of  the  11  states  which  implies  a 
locking  of  S  to  the  internuclear  axis.  This  effect  is  high¬ 
lighted  by  applying  a  frame  transformation  of  the 
IT M ,p\£,f,S ,/)  channel  basis  such  that  11  becomes 
represented  as  a  ‘  ‘good’  ’  quantum  number.  However,  in 
this  new  basis  /  and  consequently  ('  are  no  longer  con¬ 
served.  Fortunately  the  splitting  of  the  degeneracy  of 
the  a  3£u  state  is  small,  and  hence  /  and  £  still  remain 
good  approximate  quantum  numbers.  This  allows  us  to 
block  the  Hamiltonian  for  a  given  F  and  p  into  subspaces 


of  common  f  and  £  values.  The  atomic  stretched  state 
angular  momenta,  /a  =/b  =  2  and  md  =  mb  =  2,  can  only 
couple  to  an  /=  4  state,  which  then  couples  with  the 
£  =  0  s  -wave  to  give  a  total  angular  momentum  F  =  4. 
This  is  the  only  F  value  for  which  there  is  a  stretched 
state  ,v  -wave.  The  spin-spin  dipole  interaction  can  only 
change  £  by  two  units,  and  thus  upon  examining  the 
F  =  4,  p  =  +  1  Hamiltonian,  we  find  that  the  (f=  4, 
£=0 )  subspace  can  only  couple  directly  to  the  (f=  3, 
£  =  2)  and  (f=  2,  £=  2)  subspaces. 

Our  calculations  are  carried  out  with  the  full  compli¬ 
ment  of  accessible  channels  associated  with  a  given 
F ,M,p  and  our  rates  are  obtained  with  a  sufficient  sum¬ 
mation  over  F,M,p  to  insure  convergence.  However,  we 
find  that  at  the  temperatures  relevant  to  BEC  (T  <  100 
p,K)  only  the  single  set  of  F  =  4,  p  =  +  1  solutions, 
which  involve  the  close  coupling  of  20  channels  in  Eq. 
(3),  contribute  significantly  to  the  stretched  state  spin 
relaxation.  This  is  because  only  the  incident  ,v  -wave  con¬ 
tributes  to  spin  relaxation,  since  the  contributions  from 
incident  channels  with  /  >  2  are  strongly  suppressed  at 
these  temperatures  due  to  quantum  threshold  effects. 
Furthermore,  of  the  20  channels  contributing  to  Eq.  (3) 
for  F  =  4,  p  =  +  1,  and  coupled  by  the  20  X  20  interac¬ 
tion  matrix  Uyy(R),  only  the  five  diabatic  channels 
labelled  y  =  1-5  in  Table  1  and  consisting  of  the  three 
subspaces  described  above,  play  any  significant  role. 

Figure  la  shows  the  diagonal  interaction  potentials 
Uy,y(R)  for  these  five  channels.  Because  of  complicated 
curve  crossings  and  strong  interactions  in  the  y  basis, 
more  insight  comes  when  we  examine  the  five  a  =  1 ,5 
“adiabatic”  potentials  Va(R)  shown  in  Figs,  la  and  lb. 
These  are  obtained  by  diagonalizing  the  5X5  interac¬ 
tion  potential  Ury(R)  at  each  R.  The  five  potentials 
never  cross  and  are  labeled  in  order  of  increasing  energy. 
Each  adiabatic  channel  a  correlates  asymptotically  with 
the  corresponding  y  state  in  Table  1 .  At  short  R  each  a 
potential  corresponds  to  a  very  good  approximation  with 
either  the  A  'X/  or  a32+u  potentials.  Thus,  at  short  R , 
and  out  to  where  the  exchange  interaction  term  Cexc  in 
Eq.  (2)  remains  dominant,  the  Va=i(R)  potential  essen¬ 
tially  mimics  the  pure  A'S+g  potential,  and  the  other  four 
Va(R )  are  basically  pure  a 3X+U  potentials. 


Table  1.  Significant  F  =  4,  p  =  +  1  channel  states  for  spin-depolarization  rates  near 
threshold 


Diabatic  (asymptotic)  basis 

Adiabatic  Basis 

y 

(  f 

/a 

h 

Asymptotic  energy  (mK) 

a 

Short  R  label 

t 

2  2 

1 

t 

-656.022 

i 

‘s+ 

8 

2 

2  2 

2 

t 

-328.011 

2 

32+„ 

3 

2  3 

2 

i 

-328.011 

3 

32+u 

4 

0  4 

2 

2 

0.000 

4 

3r„ 

5 

2  2 

2 

2 

0.000 

5 

32+u 
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Fig.  la.  Diagonal  potentials  for  the  five  most  significant  F  =  4, 
p  =  +  1  channels.  The  dashed  curves  show  diagonal  elements  of  the 
y  =  1-5  case(e)  diabatic  interaction  matrix  U%y(R).  The  solid  curves 
show  the  adiabatic  potentials  Va(R)  for  a=  1-5  in  Table  1.  On  the 
scale  of  this  figure  the  two  y  =  3,4  and  the  four  a  =  2, 3,4,5  channels 
essentially  track  the  a3X+u  potential.  The  y  =  1,2,5  potentials  are  ad¬ 
mixtures  of  Vo,o  and  V\tn  and  are  strongly  coupled  by  off-diagonal 
terms  proportional  to  the  exchange  term  in  Eq.  (2).  The  adiabatic 
a  =  1  channel  tracks  the  X!X+g  potential. 


Fig.  lb.  Asymptotic  correlations  of  five  adiabatic  potentials  a  =  1-5 
in  Table  1. 

Solving  the  set  of  five  close-coupled  equations  and 
obtaining  the  5X5  scattering  matrix  S (y,y')  =  S(a,a') 
for  this  abbreviated  set  of  channels  is  sufficient  to  quan¬ 
titatively  reproduce  the  elaborate  multichannel  calcula¬ 
tions  to  within  a  few  percent  at  all  energies  below 
100  piC  What  has  been  defined  [8]  as  the  “event”  rate 
constant  Advent  is  simply  the  sum  over  all  inelastic  events 
experienced  by  the  incident  stretched  state  channel 
y  =  a  =  4, 


Keve =  (^  )  [ISl/  +  ISl 4f  +  IS3/  +  IS5.4l2]  (4) 

where  /j.  is  the  reduced  mass  of  the  Rb2  dimer,  and 
k  =  V(2  iie/fi2)  is  the  wave  number  of  the  a  =  4  channel 
incident  with  relative  kinetic  energy  e.  At  threshold  the 
first  three  terms  I5„,4I2  a=  1,2,3  vary  as  k  and  their 
contribution  to  K„cm  approaches  a  constant  at  low  ener¬ 
gies.  These  inelastic  elements  measure  the  coupling  to 
the  exothermic  channels  and  invariably  cause  loss  of  two 
Rb  atoms  from  the  trap.  The  fourth  term  produces  dis¬ 
orientation  of  the  stretched-state  atoms  which  ultimately 
leads  to  loss  from  the  trap  as  well.  However,  as  this 
element  vanishes  as  15s, 4I2  ^k2,  its  contribution  is  negligi¬ 
ble  as  k—>  0  and  can  be  neglected. 

The  close  coupled  equations  can  be  numerically 
solved  using  either  the  a  or  y  basis.  Since  these  basis 
sets  are  asymptotically  equivalent ,  they  lead  to  exactly 
the  same  5-matrix.  In  actual  practice  the  y  basis  is 
vastly  more  convenient  in  solving  the  close-coupling, 
while  the  a  basis  is  much  more  useful  in  gaining  physi¬ 
cal  insight  from  the  results.  For  example,  examining  the 
non-adiabatic  coupling  [17]  we  find  that  only  channels 
a  =  1  and  a  =  2  are  strongly  coupled,  and  this  occurs 
over  a  very  limited  region  R  ~  20  aa  to  25  aa  (see  be¬ 
low).  All  the  remaining  couplings  are  weak  and  pertur¬ 
bative.  This  latter  feature  will  play  an  important  role  in 
assessing  the  sensitivity  of  the  rates  to  the  singlet  poten¬ 
tial  in  Sec.  4. 

3.  Sensitivity  Analysis  of  the  Triplet 
Scattering  Length 

All  our  rates  are  calculated  using  molecular  potentials 
obtained  from  ab  initio  MCSCF  codes  employing  the 
highly  accurate  ab  initio  pseudopotentials  of  Krauss  and 
Stevens  [10]  for  R  <  20  a0,  and  joined  on  to  the  well-an¬ 
alyzed  long  range  dispersion  potentials  for  the  diatoms 
[11]  for  R  >  20  a0  Since  the  calculation  of  ultracold 
collision  rates  requires  extraordinary  accuracy  of  the 
molecular  potentials,  these  potentials  can  only  serve  as 
excellent  initial  estimates  of  the  short  range  portions  of 
the  Vs,n(R)  in  Eq-  (2).  In  particular,  they  are  not  accu¬ 
rate  enough  to  confidently  calculate  threshold  properties 
such  as  the  scattering  length.  These  parameters  are  de¬ 
termined  by  an  integration  of  Vs,n(R)  over  the  entire 
range  of  R . 

Since  we  have  confidence  in  the  long  range  parame¬ 
ters  in  Eq.  (2),  we  only  vary  the  short  range  potential  in 
order  to  assess  the  sensitivity  of  the  scattering  calcula¬ 
tions  to  these  potentials.  We  have  added  an  adjustable 
harmonic-like  short  range  term  to  the  fl32u+  and  Al2g+ 
potentials  as  follows. 
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V\0(R )  =  V0,o(R )  +  C0(R  -  Re 0)2  R<RS „= 7.89032  a0 

(5a) 

Vsri  ,a(R )  =  Vi  ,n(R )  +  Ct(R~  Rel)2  /?<7?el=  1 1 .6683  a0 

(5b) 

where  Eq.  (5)  is  only  applied  at  distances  smaller  than 
the  indicated  equilibrium  internuclear  distances  for  the 
attractive  singlet  and  triplet  potentials  respectively.  Note 
that  the  functional  form  in  Eq.  (4)  is  chosen  arbitrarily 
and  has  no  physical  significance;  other  short-range 
forms  could  be  used  with  equal  effect.  In  particular  we 
want  to  show  that  by  choosing  an  appropriate  value  for 
Ci  in  Eq.  (4b)  we  insure  that  the  associated  Vi,n(R) 
potentials  give  a  scattering  length  of  1 10  ao  as  suggested 
by  the  joint  theoretical-experimental  analysis  of  Boesten 
et  al .  [16]. 

Our  initial  fit  to  the  ab  initio  calculations  [10]  for  the 
a  ’’X  potential  happened  to  support  38  vibrational 
levels.  The  last  level  v  =  37  has  a  binding  energy 
of  -  1.45618  X  10“8  au  (1  au  =  e2/a„  =  4.359748  X 
10~18  J,  where  e  is  the  electron  charge),  and  a  positive 
scattering  length  of  21.6  aa.  This  corresponds  to  the 
point  at  Ci  =  0  in  Fig.  2  where  we  show  the  variation  in 
scattering  length  as  we  systematically  vary  Ci  in  Eq. 
(5b).  The  figure  shows  that  the  scattering  length  A,  is 
extremely  sensitive  to  the  short-range  portion  of  the 
potential.  The  scattering  length  is  defined  by  the  .5-wave 
threshold  behaviour  of  the  elastic  scattering  phase  shift 
£i  — >  —  kA  i  in  the  limit  that  the  asymptotic  wavenumber 


Fig.  2.  The  variation  of  s7Rb2  scattering  length  Ai  for  the  entrance 
s-wave  channel  incident  on  the  a3£+u  potential,  as  a  function  of  the 
short  range  Ci  parameter,  extracted  from  the  wavefunction  at 
e/kB  =  0.1  nK  (full  circles).  The  solid  line  through  the  points  gives  the 
fit  of  Eq.  (6).  We  also  plot  the  threshold  values  v(0)  of  the  bound  state 
phase  as  we  vary  Ci  (full  squares).  The  modular  ir  values  of 
v(e„)  =  rrn  identify  the  bound  state  eigenvalues.  As  the  last  bound 
state  e„  eigenvalue  approaches  zero  the  predicted  scattering  length 
passes  from  a  large  positive  value  to  a  large  negative  value  as  the  level 
"pops"  out  of  the  potential.  This  occurs  at  about  Ci  =  —  2.6  X  10 
au/flo2  for  87Rb2  where  v(0)  — >  38ir. 


k  goes  to  zero  (frk~l 2p  =  e  with  e  the  collision  energy). 
It  is  well  known  [20]  that  the  actual  value  and  the  sign 
of  At  is  critically  dependent  on  the  position  of  the  last 
bound  state  that  can  be  supported  by  a  given  potential. 
This  in  turn  is  related  to  what  we  like  to  call  [17,21]  the 
bound  state  phase  v(k),  which  is  defined  for  negative 
energies  e  =  —  h1  k2  12  p,  which  lie  below  the  threshold  at 
e  =  0,  and  where  k  in  turn  is  defined  as  a  continuous 
positive  real  variable.  The  modular-TT  value  of  the  bound 
state  phase  v(k„)  =  «tt  locates  the  position  of  the  nth 
vibrational  eigenvalue  [21,25],  such  that  e„  =  —fi2Kn2l 
2/jl.  Actually  the  deviation  of  the  threshold  value  of  v(0) 
from  modular-Tr  is  a  useful  measure  of  the  last  bound 
state  position.  This  quantity  plays  a  prominent  role  in 
many  descriptions  of  threshold  behaviour  (see  Stwalley 
[22],  LeRoy  and  Bernstein  [23]  and  the  Eindhoven 
group  [16]),  where  it  has  been  denoted  as  Vd  and  is 
sometimes  called  the  effective  vibrational  quantum 
number  at  the  dissociation  limit.  For  our  initial  fit  with 
Ci  =  0  we  found  v(0)  =  37.699541tt.  This  quantity  in¬ 
creases  or  deceases  monotonically  as  we  systematically 
make  the  potential  more  or  less  attractive  by  varying  Ci 
in  Eq.  (4b).  At  9.75  a0,  the  zero  energy  turning  point  of 
the  a  'X+u  potential  for  Ci  =  0,  a  value  of  Ci  = 
±5  X  10-5  au/a02  produces  a  ±  50  cm-1  change  in  the 
potential.  For  comparison,  the  triplet  state  potential  is 
204  cm-1  deep  at  its  potential  minimum  Re i. 

Figure  2  shows  that  the  scattering  length  as  a  function 
of  the  Ci  shift  parameter  passes  from  plus  to  minus 
infinity  as  the  last  bound  state  is  pushed  out  of  the 
potential.  The  position  of  the  singularity  is  easily  lo¬ 
cated  by  examining  the  threshold  behaviour  of  the  bound 
state  phase  v(0).  We  see  that  this  quantity  approaches 
38tt  just  as  Ci  approaches  —  2.6  X  10-5  au la02  and  Ai 
passes  through  infinity.  In  fact  a  nice  analytic  relation¬ 
ship  exists  between  the  scattering  length  A  and  v(0). 


A=-f 


cot 


+  cot  v(0) 


(6) 


where  e  =  — 


and  the  derivative  is  defined  by 
k.  The  parameter  s  in  cot(TT/ 


h2K2 

2p 

v(K)  =  v(  0)  +  !^ 

(s  —  2))  is  defined  by  the  leading  asymptotic  power  law 
—  CJR '  for  the  potential.  Both  the  a32+u  and  the  A  'S+g 
potential  have  as  the  leading  term  CeR  6  and  we  have 
cot(TT/(s  —  2))  =  1.  (Instead  of  the  pme  a3X+u  potential 
we  actually  prefer  to  use  the  adiabatic  potential  desig¬ 
nated  as  a  =  4  in  Table  1,  in  which  case  the  lower  order 
a  R  3  terms  in  Eq.  (2b)  and  (2c)  are  rigorously  removed 
by  the  diagonalization  of  the  interaction  matrix  and  do 
not  contribute  to  the  threshold  behaviour  of  this  s -wave 
channel).  If  we  evaluate  v(e)  at  an  eigenvalue  e  =  e„  then 
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v  (0)  ~  hjt  — —  Kn  =  htt  +  8„.  In  the  special  case,  as 

O  K 

e„  -A  —  0  and  the  last  bound  state  lies  just  below  the 

dissociation  limit,  such  that  8„  ~  —  -r—  k„  «  1  and 

c)k 

tanS„  -A  S„,  we  obtain  the  usual  perturbative  expression 
(Ref.  [20]  p.  48)  for  the  scattering  length 


A 


dv 

c)k 


dv  J_  J_ 

3  k  8„  k„  ' 


(7) 


Although  the  conventional  derivation  of  this  expression 
is  limited  to  a  potential  with  a  single  bound  state,  we 
find  perfect  agreement  with  this  limiting  behavior  even 
for  wells  supporting  many  bound  states. 

Actually,  if  we  have  a  situation  where  v(0)  is  not  quite 
ready  to  support  the  n th  bound  state,  i.e.,  v (0)  =  m r  —  8„ 

such  that  8„  =  —  -r—  (k„)  «  1,  then  we  can  visualize  a 
3k 

“pseudobound”  state  lying  just  above  the  dissociation 
limit  with  an  “eigenvalue”  en  =  +  A1  K;J2j±.  In  the  limit 
where  tan  v(0)  -4  —  8„  we  obtain  an  expression  which 

complements  Eq.  (7),  A  — >  —  i  ,  and  predicts  a  nega- 

K„ 

five  scattering  length  whenever  a  pseudobound  state  lies 
just  above  the  dissociation  limit.  This  behavior  is  well 
substantiated,  and  quantitatively  confirmed,  by  the  re¬ 
sults  in  Fig.  2. 

Numerical  studies  show  that  the  bound  state  phase 

dv 

does  indeed  vary  as  v(k  )  =  v(0)  +  —  k  near  threshold, 

dv 

and  furthermore  the  quantity  —  required  in  Eq.  (6)  is 

basically  an  asymptotic  property  that  only  depends  on 
the  long-range  potential,  and  is  totally  insensitive  to 
variations  in  the  short  range  potential,  such  as  those 
introduced  by  the  shift  parameter  C\.  From  these  data 
dv 

we  estimate  that  —  ~  —78  a0.  Using  this  estimate  we 
3k 

have  plotted  Eq.  (6)  as  the  solid  curve  in  Fig.  2  and  find 
perfect  agreement  over  the  entire  (modular  tt)  range  of 
v(0)  with  the  calculated  points.  Note  that  this  expression 
also  predicts  the  exact  locations  of  the  zeros  in  the 
scattering  length,  which  are  always  located  at 
v  (0)  =  n  7T  +  0.75tt.  It  is  interesting  to  note  that  over  the 
modular  tt  range  the  scattering  length  is  predicted  to  be 
positive  for  3/4  of  the  range,  and  negative  for  only  1/4. 
This  means  that,  if  we  know  nothing  about  the  short 
range  potential,  we  can  at  least  predict  there  is  a  3:1 
probability  that  the  scattering  length  will  be  positive! 
The  functional  form  of  Eq.  (6)  was  also  confirmed  for 
an  asymptotic  R  8  by  setting  C(,  =  0  in  Eq.  (2),  such  that 
cot(rr/(s  —  2))  =  cot(Tr/6)  =  V3. 


If  we  choose  C\  =  3.128  X  10~5  au/a02  we  obtain  a 
scattering  length  AJ87)  =  +  109.1  a0  and  a  threshold 
value  of  Vg7(0)  — >  37.3761rr,  which  predicts  38  bound 
vibrational  levels  with  a  last  bound  state  at  E 37  = 
—  2.947925  X  10-9  au.  This  choice  is  made  to  conform 
to  the  scattering  length  obtained  by  the  Texas/Eindhoven 
group  [15,16].  Further  confidence  is  obtained  from  Fig. 
3  where  we  compare  the  scattering  lengths  for  the  S7Rb 
and  the  85Rb  isotopes.  In  calculating  the  two  scattering 
lengths,  the  only  difference  is  the  mass  of  the  two  iso¬ 
topes.  We  see  that  at  the  same  value  of  C\  the  calculated 
scattering  length  Ai(85)  =  —  309. 1  a0  for  85Rb  is  in  good 
agreement  with  [15,16].  In  addition,  the  threshold  value 
V85(0)  =  36.9367tt  predicts  37  vibrational  levels  with  a 
last  bound  state  at  E36  =  —  3.305245  X  10-8  au. 


Ct  (10  5  au  aQ  2) 

Fig.  3.  Scattering  lengths  as  in  Fig.  2,  for  the  s7Rb(solid  circles)  and 
85Rb(open  circles)  isotopes.  These  are  evaluated  for  the  same  a32+„ 
potential,  using  the  appropriate  mass  for  each  isotope.  For  all  future 
calculations  we  use  C i  =  3.13  X  10  au/<v,  indicated  by  the  vertical 
arrow.  This  choice  yields  scattering  lengths  Ai(87)  =  109.1  a0  and 
Ai(85)  =  —  309.3  «0,  as  prescribed  by  the  analysis  given  by  the  Texas/ 
Eindhoven  group  [15.16], 


One  final  confirmation  of  the  validity  of  the  a 
potential  we  are  using  is  obtained  by  examining  the 
d- wave  shape  resonance  structure  defined  by  the  adia¬ 
batic  channel  a  =  5  in  Table  1 .  This  corresponds  to  an 
incident  channel  which  correlates  with  the  triplet  state  at 
short  distance,  and  correlates  with  th e/a  =/b  =  2  atomic 
states  at  large  distance  entering  with  an  £  =  2  partial 
wave.  This  gives  rise  to  the  centrifugal  barrier  shown  in 
Fig.  4,  with  a  barrier  height  of  420  |xK.  The  radial 
wavcfu lie tio u s  /„=4( R )  iind  f^AR  )  are  shown  in  Fig.  5  for 
four  incident  kinetic  energies:  100  p.K,  200  p,K,  350 
p,K,  and  500  p,K.  These  f„(R)  are  the  single-channel, 
energy  normalized  elastic  scattering  wavefunction  asso¬ 
ciated  with  the  Va(R )  adiabatic  potential.  Boesten  et  al. 
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Fig.  4.  The  ( =  2  centrifugal  barrier  for  the  a  —  5  channel,  which 
asymptotically  correlates  with  the  a3£+n  potential.  The  four  energies 
indicated  will  be  used  in  Fig.  5. 


Fig.  5.  Energy  normalized  continuum  wavefunctions  for  the  adiabatic 
a  =  5  channel  for  energies  e  =  100  p.K,  200  p,K,  350  p,K,  and  500  jx K 
(dashed  lines)  and  for  the  adiabatic  a  — 4  channel  for  energy  e  =  350 
jx K  (solid  line).  The  adiabatic  a  =  5  channel  exhibits  d- wave  shape 
resonance  structure  with  maximum  amplitude  enhancement  near 
Emax  =  350  p.K. 


[16]  have  concluded  that  there  is  a  shape  resonance 
enhancement  of  the  photoassociation  from  channel 
a  =  5  relative  to  the  a  =  4  channel  with  incident  ('  =  0. 
The  solid  curve  shows  the  a  =  4  wavefunction  for  350 
|xK.  The  a  =  5  channel  amplitude  increases  with  energy 
up  to  e  ~  350  p.K,  after  which  it  decreases.  One  would 
infer  from  these  plots  that  the  d -wave  shape  resonance 
is  very  broad  with  a  ‘  ‘width’  ’  extending  from  200  |xK  to 
500  |xK.  This  behaviour  is  consistent  with  the  analysis 
given  in  reference  [16]. 


4.  Sensitivity  of  Loss  Rates  to  the  Singlet 
Potential 

As  we  systematically  varied  the  short-range  A'X+g 
potential  by  varying  Co  in  Eq.  (5)  we  find  a  surprisingly 
large,  and  seemingly  erratic  variation  in  the  stretched 
state  loss  rate  coefficient.  The  spin-relaxation  rates  are 
determined  by  the  very  small  5 (a, 4)  matrix  elements, 
which  typically  have  a  magnitude  of  10-3.  We  found  that 
the  variation  was  intimately  associated  with  the  varia¬ 
tion  of  the  large  I  .S' (a  =  1  ,a  =  2)l2  element  shown  in  Fig. 
6.  As  seen  from  Fig.  1  and  Table  1,  this  5-matrix  ele¬ 
ment  measures  the  probability  for  the  /a,/b  = 
1,2  — >/a,/b=  1,1  transition.  The  coupling  which  deter¬ 
mines  this  transition  probability  is  the  short-range  ex¬ 
change  potential,  which  causes  the  strongest  and  only 
non-perturbative  inelastic  event  associated  with  the 
channels  in  Table  1.  Since  5(1,2)  is  evaluated  at  a  total 
energy  determined  by  the  a  =  4  entrance  channel, 
namely  0. 1  |xK  above  the  a  =  4  channel  threshold,  the 
asymptotic  kinetic  energies  in  channels  a  =  1  and  a  =  2 
are  656  mK  and  328  mK,  respectively  (see  Table  1). 
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Fig.  6.  Variation  of  IS12I2  (for  a  total  energy  e/kB  =  0.1  p.K  incident  in 
channel  4)  as  a  function  of  the  short-range  Co  parameter  (dashed  line) 
and  variation  of  v(0)  as  a  function  of  the  same  parameter  (dotted  line). 
At  Co  =  8.0  X  10~5  au/a02  the  quantity  v(0)/ir  =  121  and  the 
potential  supports  exactly  122  vibrational  states.  At  this  point  the 
singlet  scattering  length  passes  through  infinity.  However,  IS12I2  varies 
smoothly  through  this  region,  because  channels  1  and  2  have  high 
asymptotic  kinetic  energy,  well  away  from  threshold. 


Before  we  can  understand  the  strong  influence  of  the 
ArlX+g  potential  on  the  perturbative  5 (a ,4)  elements  we 
must  first  examine  the  profound  effect  of  this  potential 
on  the  behaviour  of  5(1,2).  The  change  to  the  short 
range  potential  is  sufficient  to  make  small  displacements 
in  the  nodes  of  fa=l(R)  in  the  coupling  region  which  is 
important  for  determining  5(1,2),  but  it  is  not  immedi- 
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ately  obvious  why  this  leads  to  such  large  variations  in 
the  above-threshold  5(1,2).  Fig.  6  shows  I5(l,2)l2  varies 
from  a  minimum  of  0  at  Co  ~  —  9.4  X  10-5  au /a02  to  a 
maximum  of  0.834  at  Co  ~  —  16.4  X  10-5  au/a02. 

This  variation  is  best  understood  by  considering  the 
matrix  element  of  the  radial  coupling  operator 

Q]  2(R )  ^77  [17]  between  the  adiabatic  states  a=  1,2  in 
oK 

Table  1  and  Fig.  7.  In  the  first  order  distorted  wave 
approximation  (Ref.  [20]  p.  349,  and  Refs.  [26  and  27]) 
the  5(1,2)  matrix  element  is  proportional  to  the  integral, 

5(1,2)  °c  JdRfi(R)Ql2(R)  ~^f2(R), 

where  Qn(R)  is  determined  by  the  R -variation  of  the 
orthogonal  5X5  matrix  Mn(R )  which  diagonalizes  the 
diabatic  interaction  matrix,  U(R)  =  M(R)V(R)M(R), 
such  that 

Qn(R)  =  -=i.5  MUj(R)  dM2j(R )/dR . 

Not  surprisingly  we  find  the  coupling  is  highly  localized 
in  the  vicinity  of  R  ~  22  (see  Fig.  7)  where  the  spin-ex- 
change  splitting  between  the  singlet  and  triplet  potential 
in  Eq.  (2)  becomes  comparable  to  the  hypeifine  splitting 
of  the  atoms.  Although  Q i2  causes  strong  nonadiabatic 
mixing  between  channels  1  and  2,  the  distorted  wave 
approximation  above  is  suitable  for  the  qualitative  argu¬ 
ment  we  make  below,  even  though  it  is  not  suitable  for 
quantitative  calculations. 


Fig.  7.  The  adiabatic  potentials  for  channels  a  =  1  and  a  =  2  together 
with  the  non-adiabatic  coupling  operator  (in  arbitrary  units)  between 
these  channels.  The  latter  is  well  localized  in  the  region  of  18  a0  < 
R<  26  do,  where  the  exchange  splitting  between  the  -CSV  and  a32+u 
potentials  becomes  comparable  to  the  hyperfine  splitting. 


Changing  the  singlet  potential  has  a  negligible  effect 
on  the  adiabatic  /:(R ),  which  is  shown  by  the  solid  curve 
in  Fig.  8  and  primarily  portrays  a  pure  triplet  state,  at 
least  up  to  the  vicinity  of  strong  coupling  near  22  a„. 
The  two  different  short  range  singlet  potentials,  associ¬ 
ated  with  the  indicated  minimum  and  maximum 
15 (1 ,2)l2  in  Fig.  6,  are  used  to  obtain  the  two  different  /] 
functions  which  we  designate  as /i!min  and/i  max  in  Fig.  8. 
Note  the  structure  of  these  functions  in  the  vicinity  of 
the  peak  of  the  non-adiabatic  coupling  operator  Ql2.  The 
Amin  function  has  an  almost  perfect  overlap  with  the  /, 
function.  Since  the  complete  non-adiabatic  operator  is 
equal  to  Qi2  times  the  radial  derivative  d/dR  this  perfect 
overlap  implies  a  very  poor  overlap  between  /]  and 
dfi.mmiR  )/dR ,  and  therefore  implies  that  the  distorted 
wave  integral  should  be  quite  small.  In  fact,  for  this  case 
our  exact  close-coupling  results  yield  15  ( 1 ,2)l2  = 
0.000017  (presumably  by  varying  Co  slightly  we  could 
have  found  a  perfect  cancellation  with  15(1 ,2)l2  =  0).  The 
second  function  /2jmax  in  Fig.  8  is  phase  shifted  with 
respect  to/,,  implying  improved  overlap  between _/]  and 
’Ifi.mixJR  )/dR  and  a  larger  distorted  wave  integral.  For 
this  case  our  close-coupling  predicts  the  maximum 
15(1, 2)l2  =  0.83. 


Fig.  8.  The  adiabatic  /2  wavefunction  (solid  curve)  and  the  adiabatic 
wave  functions  f\  imax  (dotted  curve)  and/i,min  (dashed  curve).  The  latter 
were  calculated  using  the  singlet  shift  parameters  Co=  —  16.4  X  10-5 
au /aQ2  and  Co  =  —  9.80  X  10-5  au /a02,  which  produce  the  maximum 
and  the  minimum  inelastic  LS12I2  elements  in  Fig.  7  respectively.  The 
<2 12  operator  (in  arbitrary  units)  is  also  shown,  and  locates  the  region 
of  strong  non-adiabatic  coupling  between  channels  1  and  2. 

Our  initial  fit  to  the  ab  initio  calculations  [4]  for  the 
*X+g  potential,  combined  with  the  3X+U  potential  with  a 
scattering  length  Ai=  109  ao  prescribed  by  Ref.  [16], 
just  happened  to  yield  a  value  of  15 ( 1, 2)l2  =  0.0731. 
When  the  L2+g  potential  was  systematically  modified  at 
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short  distances  by  increasing  Co,  with  nodes  being  sys¬ 
tematically  pushed  to  larger  distances.  Fig.  9  shows  that 
event  tracks  the  associated  values  of  15(1, 2)P  with 
varying  by  roughly  a  factor  of  4.  The  figure  also  shows 
that  several  of  the  individual  components  contributing  to 
the  spin-relaxation  show  similar  qualitative  variation 
with  differences  in  detail. 


Fig.  9.  IS,2I:  and  Kcven,  versus  the  short-range  Co  parameter  for  the 
ground  X'£+g  potential.  The  figure  shows  the  strong  correlation  be¬ 
tween  the  two  quantities.  The  correlation  with  specific  components 
^fcvent(22,22  — >/ama,  /bmb)  are  also  shown:  (a)  /ama,  /onto  -11.11; 

(b )/ama,/bmb  =  22,10;  (c)/ama,/bmb  =  21,1 1;  (d) /ama,/bmb  =  22,1 1.  All 
rate  coefficients  have  been  divided  by  2  X  10  15  c  m  Vs  in  order  to 
place  them  on  the  same  scale  as  the  dimensionless  ISi2l2. 

Figure  10  shows  the  behavior  of  Kewnt  with  respect  to 
the  variation  in  I5(l,2)l2.  We  show  the  dependence  of  the 
stretched  state  rate  coefficients  on  the  A'X+g  potential 
for  three  different  situations.  The  SS-only  curve  gives 
Kw cm  when  we  only  include  the  usual  spin-spin  splitting 
which  varies  as  3a2/2R3.  The  SO-only  curve  shows  the 
very  small  rate  coefficient,  almost  independent  of 
5(1,2),  that  results  if  the  SS  splitting  is  removed  and 
only  the  short-range  second-order  spin-orbit  splitting  is 
included.  In  Sec.  5  we  will  show  that  the  SO-only  con¬ 
tribution  primarily  occurs  at  short  distance,  R  <  20  a0, 
where  significant  molecular  interactions  can  occur. 
Consequently,  the  dependence  of  the  SO-only  rate  coef¬ 
ficient  on  15(1, 2)l2  is  minuscule.  Finally  the  SS-plus-SO 
curve  shows  the  rate  coefficient  when  both  SS  and  SO 
interactions  are  included  in  the  close-coupling  calcula¬ 
tion.  This  curve  demonstrates  up  to  a  factor  of  two 
reduction  in  the  rate  coefficient  when  compared  to  the 
SS -curve 

The  rate  coefficients  in  Fig.  10  are  double  valued  as 
a  function  of  15 ( 1 ,2)l2  except  at  the  extremes  of  the  range. 
This  property  can  be  traced  to  Fig.  9,  where  we  have 


Fig.  10.  Rate  coefficient  Severn  versus  IS12I2  as  Co  is  varied  for  the 
potential.  The  variation  of  C0  is  sufficient  to  “pop”  two  bound 
states  out  of  the  X'X+g  potential.  Three  such  curves  are  shown.  The 
curve  labeled  “SS  only”  includes  the  spin-spin  (SS)  interaction  only. 
The  lowest  curve  labeled  “SO  only”  includes  the  second-order  spin- 
orbit  (SO)  interaction  only.  The  middle  curve  labeled  “SS  plus  SO” 
includes  both  the  SS  and  SO  interactions. 

followed  15 (1 ,2)l2  over  a  range  of  Co  that  changes  v(0)/tt 
by  unity  and  ‘  ‘pushes’  ’  one  bound  state  out  of  the  X  l2+„ 
potential  well.  We  find  two  identical  values  of  15(1 ,2)l2 
in  this  interval.  The  second  value  occurs  when  the  f\ 
function  is  shifted  by  approximately  half  a  deBroglie 
wavelength  in  the  peak  region  of  Q 12,  thereby  maintain¬ 
ing  the  same  value  of  the  distorted  wave  integral  for 
15 (1 ,2)l2  (see  Fig.  8).  This  changes  the  sign  of  5(1,2)  but 
results  in  the  same  15(1 ,2)l2.  However,  this  change  in 
sign  results  in  slightly  different  interference  effects  in 
the  evaluation  of  5 (a ,4)  and  thus  at  distances  be¬ 
yond  R  ~  26  a0,  as  we  will  discuss  below.  Actually,  in 
Fig.  10  we  varied  Co  enough  to  “eject”  two  bound  states 
from  the  potential,  and,  for  a  given  resultant  15(1, 2)l2  the 
rate  coefficients  in  Fig.  10  can  not  be  distinguished. 

At  the  extrema  of  Kevcnt  in  Figs.  9  or  10,  it  is  possible 
to  associate  a  single  A0  scattering  length  for  the  X  'S+g 
potential  with  the  particular  15 (1 ,2)l2  value.  This  is  not 
possible  away  from  the  extrema,  since  there  are  two 
values  that  correspond  to  the  same  Keve The  minimum 
in  /Cve nt  corresponds  to  a  scattering  length  for  the  X  l2+„ 
potential  of  A0  =  +  95  a0  and  the  maximum  corresponds 
to  a  value  of  A0  =  +  54  a0.  If  the  scattering  length  were 
measured  to  be  near  one  of  these  values,  then  the  relax¬ 
ation  rate  will  be  near  one  of  its  extreme  values. 

We  will  now  present  a  qualitative  argument  why  vary¬ 
ing  the  magnitude  of  5(1,2)  affects  the  magnitude  of  the 
spin-relaxation  rate  involving  the  5 (a ,4),  a  A  4,  matrix 
elements.  Since  the  coupling  is  weak  it  is  an  excellent 
approximation  to  represent  these  elements  as  follows; 
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5  (a  ,4)  oc  2  <  F~,a  I  Q,,a(R  )  ^  I A  )  18) 


d  x 

0  x 

_x  ~d 

_x  0 

(10) 


where  F~>a  are  defined  in  the  same  manner  as  in  Eq.  (3), 
but  now  in  the  adiabatic  channel  basis  and  with  an 
outgoing  state  in  channel  a .  Each  column  vector  in  the 
matrix  F~(R )  defines  the  radial  components  of  a  five- 
channel  close-coupled  outgoing  wavefunction  for  each 
of  the  adiabatic  channels  a  =  1,5.  The  function/!  repre¬ 
sents  the  incoming  wavefunction  of  the  a  =  4  adiabatic 
potential  of  Fig.  1 .  This  can  be  viewed  as  a  generalized, 
multichannel  version  of  the  Distorted  Wave  Approxima¬ 
tion  [20,26,27].  The  Q,,  a  matrix  elements  introduce  the 
weak  spin-dependent  coupling  between  adiabatic  chan¬ 
nels  4  and  a'. 

At  short  distances,  to  the  left  of  the  strong  coupling 
region  in  Fig.  8  the  solutions  F~(R )  are  simply  propor¬ 
tional  to  the  5X5  diagonal  matrix  of  adiabatic  refer¬ 
ence  functions  /, 


where,  especially  if  one  uses  an  adiabatic  representa¬ 
tion,  the  diagonal  elements  d  are  generally  negligible. 
For  the  strong  coupling  case,  15(1,2)1  =  12(1,2)1  can  have 
a  maximum  value  approaching  unity. 

For  the  simple  model  used  here,  we  use  x  =  1  for  the 
strong  coupling  case  and  x  =  0  for  the  weak  coupling 
case.  For  these  two  extreme  cases  the  function  F~  in  Eq. 
(9b)  takes  a  simple  form.  For  the  two  channels  a  =  1  and 
a  =  2  the  relevant  matrix  elements  of  F~  are: 


/zfe  ih\  e 
z/z/e-'^1  /z2e_,& 


x  =  1,  15(1,2)1=1  (11s) 


7ie_ifl 

0 


0 


x  =  0,  15(1,2)1  =  0  (llw) 


F~(R)  o cf(R)Q~‘(=  8a>fa(R)e~i(a  for  R  <  20  a0. 

(9a) 

Applying  generalized  MQDT  theory  [17]  and  its  associ¬ 
ated  half  collision  amplitude  [18]  to  the  outgoing  multi¬ 
channel  functions  at  distances  beyond  about  R  =  26  a0, 
the  exact  close-coupled  wavefunctions  can  be  repre¬ 
sented  rigorously  as  follows 


F~(R )  «  \f(R)  +  g(R)Y ]  (1  +  iYyle~i(  for  R>  26  a.u. 

(9b) 

where  g(R  )  =  8„„g„(R )  and  ga(R )  is  an  irregular  solu¬ 
tion  for  the  adiabatic  V„(R )  potential.  Alternatively,  this 
may  be  written  using  running  wave  reference  functions: 


^■=2 


h+  —  h  2* 


/  /z±  =  g  ±  if  — »  k  ""e 


l/2„±i(*R  +  fl 

(9c) 


where  the  real  symmetric  Y  matrix  is  related  to  the 
close-coupled  scattering  matrix  5  =  e'f2e'^=  e'-(  1  +  iY ) 
(1  —  ;Y)-1e'f.  The  only  Y  matrix  element  of  any  magni¬ 
tude  in  this  ultracold  five  channel  system  is  Y\X,  there¬ 
fore,  we  need  only  consider  the  a  =  1  and  2  channels 
and  we  can  reduce  F~  to  a  2  X  2  matrix  for  these  two 
channels.  For  the  weak  coupling  case  in  Fig.  8  even  the 
K],2  element  is  negligible  and  the  structure  of  F~(R) 
remains  diagonal  as  in  Eq.  (9a)  for  all  R .  However,  if  the 
chosen  X'2+g  potential  yields  a  large  5(1,2),  the  struc¬ 
ture  in  Eq.  (9b)  can  strongly  influence  Eq.  (8).  It  is  not 
a  bad  approximation  [17]  to  represent  a  2  X  2  Y  matrix 
as  follows. 


We  designate  these  two  cases  “s”  and  “w”  respec¬ 
tively,  for  strong  and  weak  inelastic  scattering  probabil¬ 
ity  measured  by  5(1,2).  Using  these  limits  in  Eq.  (8)  we 
find 

5(1,4)  cc-l  (h^-‘(-'\V<A\f4)  +  </z2-e-/f>IV2,4!/4> 

5(2,4)  cc-l  (h+2e-ib\V2A\f4)  +  (/zre-'^IVu!/)  (12s) 
or 


5(1,4)  cc  ft  e-ttl VlA\f4) 

5(2,4)  cc  <f2  e-ifeIV2;4!/>  (12w) 

The  radial  functions  /  and /2  are  just  the  elastic  scatter¬ 
ing  standing  waves  shown  in  Fig.  8  and  oscillate  strongly 
against  the  standing  wave  f4  defined  by  the  incident 
channel.  Thus,  we  expect  small  values  for  the  matrix 
elements  in  Eq.  (12w).  In  the  strong  coupling  regime  the 
matrix  elements  involve  the  overlap  of  /  with  pure  out¬ 
going  (or  incoming)  running  waves  with  amplitudes 
which  do  not  oscillate  with  R ,  and  we  can  easily  under¬ 
stand  why  the  matrix  elements  in  Eq.  (12s)  are  much 
larger  than  those  in  Eq.  (12w).  The  enhancement  of  the 
spin-relaxation  rate  coefficient  in  the  presence  of  strong 
final  state  interactions  is  well  demonstrated  in  Figs.  9 
and  10. 

Further  insight  into  this  final  state  effect  can  be 
gleaned  by  evaluating  K„ent  while  systematically  limiting 
the  range  of  the  spin-spin  coupling  in  Eq.  (2).  We  do  this 
by  assuming  the  fine  structure  constant  a  has  its  usual 
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constant  value  1/137  up  to  some  cutoff  distance  Rcut  and 
then  vanishes  identically  for  R  >  Rcut.  The  results  of  such 
a  study  are  shown  in  Fig.  11.  The  solid  curves  show  the 
result  using  the  X 1  X+g  potential  which  yields  the  maxi¬ 
mum  I5(l,2)l2.  In  this  case  we  expect  Eqs.  (10s)  and 
(11s)  to  prevail  at  distances  larger  than  the  peak  of  the 
Q 12  operator,  which  is  shown  by  the  dotted  curve  in  Fig. 
1 1 .  The  dashed  curves  show  the  corresponding  results 
for  the  minimum  I5(l,2)l2  case,  for  which  Eqs.  (llw) 
and  (12w)  should  be  valid  for  all  distances.  There  are 
two  sets  of  curves.  Each  set  has  one  line  with  the  SS 
interaction  only  and  one  with  our  best  estimate  of  the 
second-order  spin-orbit  terms  included.  Since  these  SO 
terms  are  short  ranged  and  their  principle  contributions 
occur  at  distances  to  the  left  of  the  <2 12  operator,  we  did 
not  apply  a  cutoff  to  these  terms. 


Fig.  11.  Kevcnt  versus  the  artificial  cutoff  parameter  Rc „,ss  of  the  spin- 
spin  interaction.  The  Vss (R )  coupling  is  set  to  zero  for  R  >  Rcm.  Kmnt 
is  shown  for  the  two  Co  values  which  yield  the  largest  (solid  lines)  and 
smallest  (dashed  lines)  inelastic  ISi2l2.  Two  curves  are  shown  for  each 
choice  of  Co.  One  corresponds  to  including  the  SS  interaction  only, 
and  the  other  corresponds  to  including  both  the  SS  and  SO  interactions 
(the  SO  coupling  in  not  cut-off).  The  arrows  show  the  values  for 

Rcu,  -> 

Both  dashed  curves  in  Fig.  11,  associated  with  the 
standing-wave  weak  coupling  solutions  in  Eq.  (llw), 
have  achieved  their  asymptotic  values  by  about  R  ~ 
40  c/,,-60  a0-  Any  possible  contributions  to  the  loss-rates 
beyond  this  region  are  quenched  by  standing  wave  oscil¬ 
lations  in  the  matrix  elements  Eq.  (12w).  On  the  other 
hand,  as  shown  by  the  behavior  of  the  pair  of  solid 
curves,  running  wave  terms  in  the  strong  field  Eq.  (12s) 
continue  to  add  to  the  rates  out  to  distance  beyond 
R  =  200  r/„  and  finally  reach  full  convergence  at  about 
R  «  300  a0. 


We  now  discuss  the  interesting  effect  that  the  total 
loss  rates  are  decreased  when  we  include  the  Vs0  terms 
as  well  as  the  V/SS  terms.  The  distorted  wave  approxima¬ 
tion  in  Eq.  (8)  implies  that  the  S (4 ,a)  S -matrix  elements 
can  be  separated  into  an  SS  and  an  SO  term.  The  loss 
rate  depends  on  the  square  of  5-matrix  elements: 

\Sss  +  Sso  I2  =  I5SSI2  +  I5S0I2  +  Sss5so*  +  Sss*Sso. 

The  first  two  terms  sum  the  individual  contributions  to 
the  rate,  whereas  the  last  terms  exhibit  the  interference 
between  them.  Figures  10  and  11  show  that  I5S0I2  is 
independent  of  the  strength  of  the  exit  channel  coupling 
measured  by  I5(l,2)l2  and  is  always  small  compared  to 
I5SSI2.  On  the  other  hand,  I.S'ssl2  does  depend  very 
strongly  on  the  strength  of  the  exit  channel  coupling. 
Obviously  the  interference  effect  is  more  significant  for 
the  case  of  weak  SS  coupling  than  for  the  case  of  strong 
SS  coupling.  In  the  former  case,  Kewm  is  decreased  by  a 
factor  of  2  when  SO  coupling  is  included,  whereas  in  the 
latter  case,  the  decrease  is  only  20  %.  Including  the  SO 
terms  causes  a  decrease  in  A'cvem,  since  the  SO  coupling 
has  an  opposite  sign  from  the  SS  coupling,  as  discussed 
in  the  next  section. 


5.  Evaluation  of  Second-Order  Spin-Orbit 
Coupling 

The  doubly  spin  polarized  atomic  states  on  the  left 
hand  side  of  Eq.  (1)  can  only  be  relaxed  by  the  weak 
coupling  between  the  two  electron  spins.  Two  terms, 
shown  in  Eq.  (2),  contribute  to  the  effective  spin-spin 
interaction  Hamiltonian:  the  direct  spin-spin  dipole  in¬ 
teraction  Vssn(R),  varying  at  long  range  as  \/R  \  and  the 
indirect  second-order  spin-orbit  interaction  Vson(R )■ 
The  latter  originates  when  the  atomic  charge  clouds 
overlap  as  a  molecule  is  formed,  and  the  interaction 
between  the  ground  state  spins  are  modified  due  to 
couplings  mediated  through  distant  excited  electronic 
states  of  the  molecule.  These  interactions  are  well 
known  in  molecular  spectroscopy  and  mimic  the  direct 
spin-spin  coupling  for  a  3S  state  [13,14],  The  reason 
these  interactions  mimic  spin-spin  coupling  is  that  they 
split  the  0  =  0  and  ±  1  components  in  a  similar  manner 
as  the  direct  spin-spin  terms.  For  heavy  species  like  Rb 
and  Cs  we  will  see  that  these  indirect  terms  can  be  much 
larger  than  the  Vssn(R)  term  at  short  distances  and 
strongly  influence  the  spin-relaxation  rate. 

We  have  calculated  ab  initio  molecular  spin-orbit  ma¬ 
trix  elements  to  obtain  estimates  of  the  second-order 
spin-orbit  correction  terms  Vson(R)  in  Eq.  (2).  For  the 
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a3Xu+  state  these  terms  are  mediated  through  distant 
electronic  states  of  'IT,  3ITU,  5IIU,  and  32u“  symmetry. 
We  calculate  the  P  -dependent  spin-orbit  matrix  ele¬ 
ments  for  these  states  using  all-electron  wavefunctions 
for  the  molecular  states  generated  by  standard  ab  initio 
methods.  The  second-order  interaction  for  the  a  ’’X„+ 
state  is  dominated  by  the  matrix  elements 
( a3Sifll//so(^)l2S+1nu,fi  )  involving  ‘IT  and  3I1U  states 
which  correlate  to  the  first  excited  2S  +  2P  atomic 
asymptotes.  In  this  case  the  second-order  coupling  term 
due  to  a  specific  2S+1nu,fl  state  takes  the  form  [14] 


V??(S,R)  =  bso(SM)Ps(R)  (13a) 


where 


bso(S,  0)  =  0, 

II 

+  1 

CO 

o 

for  5  =  0 

(13b) 

bso(S,  0)  =  2, 

cs- 

GO 

o 

? 

1+ 

II 

for  S  =  1 

(13c) 

and 


i(3-Su.n= ii//soi-5+1n„,fl= i)l 
V(2S+IU U,P)  ~  V(32U,R)  ■ 


(13d) 


Fig.  12.  Calculated  second-order  spin-orbit  coupling  parameter  of 
Eq.  (lid)  as  a  function  of  the  internuclear  separation  R.  The  circles 
show  the  contribution  from  the  first  excited  ’flu  state,  and  the  squares 
show  that  from  the  'nu  state.  Solid  points  are  for  Rb2  and  open  ones 
for  Cs2-  The  dashed  line  shows  the  difference  between  the  O  =  0  and 
lill  =  1  components  of  the  spin-spin  interaction. 


Here  HSo(R )  is  the  electronic  spin-orbit  coupling  opera¬ 
tor  and  fl  =  A  +  1 ,  where  A  is  the  projection  of  the 
electronic  orbital  angular  momentum  and  S  is  the  pro¬ 
jection  of  the  electron  spin  angular  momentum  on  the 
internuclear  axis.  Since  the  important  aspect  of  the  sec¬ 
ond-order  spin  orbit  coupling  for  dipolar  relaxation  is 
the  splitting  it  introduces  between  /2  =  0  and  ±  1  com¬ 
ponents  of  the  3Xu+  state,  we  represent  the  second-order 
spin-orbit  couplings  as  an  “effective”  spin-spin  cou¬ 
pling  term  Vn’(R)  as  follows, 

14°<>(P)  =  |[P,(P)  -  Pt,(R)] 

VKi±l(R)=-]-[Pt(R)-  Pf)(R)l  (14) 

Since  PS(P )  decays  exponentially  to  large  R ,  and  its 
magnitude  at  short  R  is  small  compared  to  uncertainties 
in  the  short  range  3XU+  potential,  we  can  ignore  the 
mean  contribution  of  these  interactions  which  equals 
[3Pt(R)  +  Po(R)]/2  and  assume  that  our  adjustment  of 
the  3Su+  potential  to  fit  the  experimental-theoretical  esti¬ 
mate  of  the  scattering  length  in  Sec.  3  actually  incorpo¬ 
rates  this  mean  spin-orbit  contribution. 

Figure  12  shows  our  calculated  Po(R)  and  !\(R )  for 
Rb2  and  Cs2.  The  b3IIu  state,  which  is  energetically  clos¬ 
est  to  the  32u+  state,  has  the  largest  coupling,  and  P\  is 


an  order  of  magnitude  larger  than  P0.  Figure  12  also 
shows  that  P0  and  Pi  have  a  much  shorter  range  than  the 
1/P3  spin-spin  term.  To  a  good  approximation  we  find 
we  can  fit  the  numerical  ab  initio  results  to  the  following 
expression, 

VsVo(P)  =  -  2VSV,(P)  =  -  Cor  (15) 

where  for  Rb2  C  =  0.001252  au,  B  =  0.975  a0_1,  and  for 
Cs2  C  =  0.02249  au,  B  =  0.830  a0~l  and  Rs  =  10  a0  for 

Q~ 

both;  here  C  is  given  in  au  (1  au  =  —  =  4.359748  X 

Clo 

10-18  J).  Note  that  Vn=o  is  negative  compared  to  the 
positive  a2/R 3  spin  dipolar  term  in  Eq.  (2b)  and  at  about 
P  ~  10.33  a„  these  two  terms  exactly  cancel  in  the  case 
of  Rb2.  The  same  is  true  for  Eq.  (2c).  Since  it  is  the 
difference  between  the  potentials  from  Eqs.  (2b)  and 
(2c)  that  gives  rise  to  spin-relaxation,  and  the  overall 
magnitude  of  these  differences  is  reduced  by  the  addi¬ 
tion  of  the  spin-orbit  interaction,  we  might  anticipate  the 
spin-dipolar  relaxation  rates  will  be  reduced  accord¬ 
ingly.  However,  we  should  emphasis  that  although  this  is 
true  for  spin  aligned  Rb,  where  the  two  contributions  are 
of  similar  magnitude,  the  opposite  effect  appears  likely 
for  Cs,  since  the  magnitude  of  the  second-order  spin-or¬ 
bit  term  for  Cs2,  although  still  negative,  is  much  larger 
than  the  spin-spin  term.  The  overall  effect  of  the  spin- 
orbit  term  in  Cs2  is  to  increase  the  spin-relaxation  rate 
relative  to  that  predicted  by  spin-spin  only. 
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To  test  the  sensitivity  of  KCYS. nt  to  the  magnitude  of  the 
spin-orbit  coupling  we  have  arbitrarily  multiplied  Eq. 
(15)  by  a  factor  of  two  and  recalculated  the  loss-rates. 
Although  it  is  impossible  to  place  unequivocal  error 
bounds  on  our  calculations  of  V7j’,  it  is  unlikely  that  our 
error  in  estimating  Vso  is  as  much  as  a  factor  of  two. 
Such  a  variation  does  give  a  reasonable  bound  to  the 
possible  effect  of  the  coupling.  In  Fig.  13  we  compare 
the  rates  for  the  SS-only  calculation  to  those  calculated 
using  Vs0  and  2  X  Vs0  added  to  Vss.  This  figure  also 
summarizes  the  status  of  our  current  confidence  in  Ktvent 
for  Rb.  Since  we  are  not  able  to  determine  the  '2+„ 
potential  with  sufficient  accuracy  to  place  any  con¬ 
straint  on  the  magnitude  of  ISi2l2,  the  rate  coefficient  is 
spanned  by  the  range  shown  by  the  curve  labeled  “SS 
plus  SO”  in  Fig.  13.  The  loss  rate  coefficient,  2Ktvent, 
lies  between  0.4  X  10'15  cm3/s  and  2.4  X  10“15  cm3/s.  If 
our  calculated  Vs0  should  be  erroneous,  we  still  expect 
the  loss  rate  coefficient  to  lie  between  0.1  X  10~15cm3/s 
and  3.0  X  10-15  cm3/s.  The  lower  bound  is  determined 
by  the  “SSplus2  X  SO’  ’  curve,  and  the  upper  bound  by 
the  “SS  only”  curve.  It  must  also  be  remembered  that 
this  assessment  is  based  on  having  a  precise  scattering 
length  for  the  fl3S+u  potential.  If  the  estimate  of  ,4  =  110 
a0  supplied  by  the  Texas/Eindhoven  should  be  in  error, 
then  the  rate  coefficient  will  be  modified  accordingly. 


Isi>2l2 


Fig.  13.  Rate  coefficient  Kcve nt  versus  IS12I2  as  Co  is  varied  for  the 
X12+g  potential.  The  upper  two  curves  are  the  same  as  in  Fig.  10.  The 
three  curves  indicate  the  uncertainty  in  KCYe nt  due  to  the  uncertainty  in 
the  potential  and  in  the  SO  interaction.  The  three  curves  show 

three  different  cases  of  the  strength  of  the  second-order  spin-orbit 
(SO)  interaction:  the  upper  curve  has  no  SO  interaction,  the  middle 
curve  uses  our  calculated  ab  initio  values,  and  the  lower  curve  uses 
twice  the  calculated  values. 


We  wish  to  make  one  final  point  before  we  conclude 
our  discussion  of  the  spin-order  effects.  Figure  12  shows 
that  Vso  is  about  an  order  of  magnitude  larger  for  Cs2 
than  for  Rb2.  In  fact,  the  contribution  to  spin-relaxation 
of  doubly  spin-polarized  Cs  atoms  from  the  second-or¬ 
der  spin-orbit  term  will  dominate  that  due  to  the  spin- 
spin  term.  Model  calculations  for  Cs  suggest  that  spin- 
relaxation  rates  for  Cs  are  insensitive  to  the  singlet 
potential.  This  agrees  with  our  result  for  Rb  shown  by 
the  nearly  horizontal  curve  labeled  “SO  only”  in  Fig. 
10.  The  spin-relaxation  rate  for  Cs  atoms  is  likely  to  be 
much  larger  than  that  for  Rb  atoms,  but  the  actual  value 
will  still  depend  on  the  a  :'X+U  potential. 

6.  Conclusions:  Assessment  of  Current 
Accuracy  of  Calculated  Rb2  Loss  Rates 

In  conclusion,  we  have  determined  the  uncertainty  in 
the  spin-relaxation  rate  coefficient  for  stretched  state 
K7Rb  atoms  associated  with  uncertainties  in  the  molecu¬ 
lar  parameters  which  control  the  magnitude  of  the  relax¬ 
ation  cross  section.  The  stretched  state  relaxation  affects 
the  lifetime  of  the  experimentally  observed  87Rb  con¬ 
densate  [1,24].  The  condensate  described  in  [1]  remains 
dilute  enough  that  its  lifetime,  ignoring  collision  with 
background  thermal  atoms,  is  determined  by  the  binary 
stretched  state  loss  rates  [24].  Ternary  rates  will  become 
dominant  in  more  dense  condensates. 

Our  calculations  use  the  a3 SJ  or  stretched  state  scat¬ 
tering  length  A 1  =  1 1 0  an  as  required  by  the  experimen¬ 
tal  photoassociation  data  for  spin-polarized  S7Rb  atoms 
[16].  Our  calculations  show  that  the  lack  of  experimen¬ 
tal  knowledge  of  the  X  'Sg  potential  provides  the  largest 
source  of  uncertainty  in  determining  the  spin-relaxation 
rate  coefficient.  The  ,v  -wave  entrance  channel  for  colli¬ 
sion  of  the /a  =  2  +/b  =  2  stretched  state  atoms  only  cou¬ 
ples  very  weakly  to  the  two  possible  exit  channels  of  the 
spin-relaxation  process,  which  produce/)  =  1  +/b  =  2  or 
/a  =  1  +/b  =  1  separated  atoms  with  increased  kinetic 
energy.  The  uncertainty  in  spin-relaxation  rate  is  associ¬ 
ated  with  the  strength  of  mixing  between  the  two  com¬ 
ponents  in  these  exit  channels.  We  varied  the  inner  wall 
of  the  X  'Xg  potential  in  order  to  determine  the  range  of 
uncertainty.  Our  calculations  show  that  the  loss  rate 
coefficient  due  to  spin-relaxation,  2Ktwent,  is  uncertain  to 
about  a  factor  of  6,  lying  in  the  range  0.4  X  10-15  cm3/s 
to  2.4  X  10-15  cm3/s.  In  a  pure  condensate  the  rate  coef¬ 
ficient  is  simply  [18].  Doubly  polarized  Rb  may 
have  the  smallest  collisional  loss  rate  coefficient  of  any 
the  alkali  species  if  the  rate  coefficient  lies  near  the 
lower  end  of  its  estimated  range. 
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We  have  provided  ab  initio  estimates  for  the  second- 
order  spin-orbit  terms  Vn’(R)  which  contribute  to  the 
effective  spin-spin  interaction.  For  the  heavier  alkali 
atoms  these  terms  have  an  important  effect  on  the  spin 
relaxation  rate.  In  87Rb,  the  VffiR )  terms  causes  a  re¬ 
duction  in  the  rate  coefficient  for  the  collisional  relax¬ 
ation  of  stretched  state  atoms  throughout  the  possible 
parameter  space  for  varying  the  X'XJ  potential.  The 
decrease  ranges  between  a  factor  of  2  %  and  20  %, 
depending  on  whether  the  rate  coefficient  lies  near  the 
lower  or  higher  end  of  its  range  of  uncertainty.  Our 
calculations  suggest  that  the  contribution  of  the  Vn°(R ) 
terms  to  the  spin-relaxation  rate  coefficient  of  stretched 
state  Cs  atoms  will  be  much  larger  than  that  from  the 
spin-spin  dipole  term. 

Our  analysis  points  out  a  critical  need  for  more  pre¬ 
cise  determinations  of  the  X  'XJ  and  a3XJ  potentials  for 
Rb2  and  Cs2  and  of  the  spin-coupling  parameters.  Exper¬ 
imental  determination  of  Vn'(R)  could  be  approached 
by  precision  spectroscopy  on  the  fine  and  hyperfine 
structure  of  the  a3X£  state.  The  quantum  chemistry  com¬ 
munity  could  provide  more  complete  and  accurate  cal¬ 
culations  of  the  second-order  spin  orbit  interactions  as  a 
function  of  R . 
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